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Tuning the Stiffness Asymmetry of Optical Tweezers via Polarization Control 
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Optical tweezers that utilize a highly focused, linearly polarized laser beam exhibit strong trap 
stiffness asymmetry, which originates from the anisotropic field distribution in the transverse plane. 
Based on the previous study of polarization-dependent focused field distribution, we explore its 
manifestation in optical trap in terms of trap stiffness asymmetry. Our results demonstrate that 
polarization control provides a versatile tuning knob for tailoring optical potential landscape even 
in case of a strongly modified focused field in the presence of dielectric spheres larger than the 
wavelength of a trap beam. 
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I. INTRODUCTION 

A tightly focused laser beam can help overcome radi¬ 
ation pressure, resulting in the stable optical trapping of 
microscopic particles in an optical potential; this princi¬ 
ple is physically realized in instruments known as “op¬ 
tical tweezer” (OTs) (TJ |2]. OTs provide a contactless 
optical clamping method and additional means of gov¬ 
erning the system dynamics (e.g., studying interaction 
dynamics in colloidal systems [aHs] and microgyroscope 
applications [II? ]). Moreover, by harnessing the extreme 
force sensitivity such methods have various applications 
in biophysics and precision control [8]. In these appli¬ 
cations, tailoring the landscape of the optical potential 
is of significant importance, and considerable efforts have 
been devoted to engineering the optical potential in order 
to realize the modeled systems desirable for manifesting 
novel characteristics of research interest (e.g., scanning¬ 
line OTs [9]). 

Recent studies on this subject in the context of highly 
focused OTs have shown that strong asymmetry in trap 
stiffnesses HDHH], which is equivalent to asymmetric trap 
potential, is primarily due to the anisotropic electro¬ 
magnetic (EM) field distribution around a trap region. 
Rohrbach et al. illustrated the trap stiffness asymmetry 
in their OT setup based on a water-immersion objec¬ 
tive lens m- Zakharian et al used the finite-difference 
time-domain (FDTD) method to investigate polariza¬ 
tion effects on the optical potential formed in air and 
in water m- Madadi et al. explored the aberration ef¬ 
fect that is caused by the refractive index mismatch be¬ 
tween the interfaces of the media involved; more specif¬ 
ically, they investigated the trap potential asymmetry 
that arises in OTs composed of an oil-immersion objec¬ 
tive lens [12]. Marag et al. elucidated the trap poten¬ 
tial anisotropy of OT by analyzing the autocorrelation 
and cross-correlation of the center-of-mass and angular 
motion of optically trapped nanotube bundle m, and 


* Electronic address: jaiminchoi@jbnu.ac.kr Fax: +82-63-270-2806 


Jinxin et al. illustrated the two-dimensional trapping- 
potential profile using the optical bottle method [T4| . 

While the abovementioned studies pd1UT2] focused on 
the characterization of trap stiffness asymmetry in OTs 
with a linearly polarized trap beam, we address the same 
subject with the aim of achieving a control parameter 
over the stiffness asymmetry by using the polarization 
state of the trapping beam. To investigate the effect 
of the polarization state of the trap laser beam on the 
trap potential asymmetry, we constructed a custom-built 
OT setup equipped with polarization control. The input 
beam is prepared with linear, elliptic, and circular po¬ 
larization, and the resulting trap potential asymmetry 
is characterized in terms of the trap stiffnesses (spring 
constants) in the transverse plane. The trap stiffnesses 
are measured for each polarization state from the motion 
of trapped particles by using the back-focal-plane inter¬ 
ferometry (BFPI) technique [15]. Theorerical estimation 
for the experimental observations is also presented based 
on the generalized Lorent-Mie theory, which is revised to 
encompass the general polarization state of a trapping 

beam nu¬ 


ll. EXPERIMENTAL SETUP AND 
PRELIMINARY MEASUREMENT 

We constructed a custom-built OT setup equipped 
with polarization control, as shown in Fig. 1(a). The 
output of the diode laser system is fed to a tapered am¬ 
plifier, resulting in 1.6 W of output power with a tuned 
wavelength (Aq) of 780 nm. The power-amplified beam 
is guided to the experimental region through a single¬ 
mode polarization-maintaining fiber, yielding the trans¬ 
mitted power of about 250 mW. The first half-wave plate 
(HWPl) and polarizing beam splitter (PBS) are used to 
adjust the input beam power, whereas the second half¬ 
wave plate (HWP2) and the quarter-wave plate (QWP) 
define the polarization state. As benchmark conditions, 
the maximum power of the incident beam and the col¬ 
limated beam size before the water-immersion objective 
lens (UPLSAPO 60XW, Olympus) were 150 mW and 4 
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FIG. 1: (a) Schematic of the experimental setup. HWP: half¬ 
wave plate. PBS: polarizing beam splitter. QWP: quarter- 
wave plate. Obj. L: a water-immersion objective lens. LI and 
L2: condenser lens and relay lens, respectively. QPD: quad¬ 
rant photodiode. LED: light-emitting diode. The dichroic 
mirror (DM) is removed during PSD measurements, (b) The 
normalized intensity distribution for a linearly polarized (LP) 
input beam at the focal plane, (c) The normalized intensity 
distribution for a circularly polarized (CP) input beam at the 
focal plane. In (b-c), the dashed white lines outline the (e“^, 
e~^) levels and the frame ticks are given units 

of Wq. 


mm (corresponding to a filling factor of 1), respectively. 
Part of that beam power was reserved for power stabiliza¬ 
tion. The wave plates were initially aligned to produce 
linearly polarized light along the horizontal direction in 
the laboratory frame:^*’^ = Eqx . A 170-/im-thick cover 
glass (N01.5H, ZEISS) and a nominal slide glass were 
sandwiched together as sidewalls in order to construct 
the sample chamber using a 200-/im-thick, double-sided 
3M tape. The polystyrene spheres (Bangs Laboratories) 
of diameter 300 ±15 nm, with the index of refraction Up 
= 1.579 at Ao = 780 nm, were diluted and dispersed in 
deuterium oxide (D 2 O) to minimize heat-induced convec¬ 
tion. Additional procedures for hydrophobic coating of 
the glass surfaces and stabilizing surfactant were made 
following the protocol provided in Ref. m- The con¬ 
denser lens (LI) collects part of the trapping beam to¬ 
gether with the scattered light. The superposed field pat¬ 
tern in the back focal plane of LI was imaged onto the 
quadrant photodiode (QPD) via the relay lens L2, which 
is known as the BFPI technique [15] . An iris at the back 
focal plane of LI was used to optimize the sensitivity of 
BFPI by adjusting the effective numerical aperture (NA) 
of the condenser lens m, and the auxiliary components 
(LED, L3, and DM) are part of the imaging setup. The 
output signals {SVx, SVy, Vsum) of the QPD were logged 
using a high-speed digital scope (5444B, PicoScope) with 
a 5-/is sampling time; each measurement was performed 
for 2 s. 

Having a diffraction-limited spot is crucial for our re¬ 
search because it ensures exclusion of other contribu¬ 


tions, e.g., aberration effects, and trap asymmetry it¬ 
self strongly depends on the spot size. We exploited 
the scattered light from the particle to produce a tight 
spot. For this diagnostic purpose, the input trap beam 
was prepared with a circular polarization, i.e.,F^o^ ^ 
Eq{x ± iy)l^/2^ and part of the scattered light from the 
test sphere being trapped was recollected by the objec¬ 
tive lens. It passed the QWP again, mostly resulting in 
the polarization that is rotated by 90°: Eo{x — iy)/\/2 
Eoy. The rotated scattered light was redirected by the 
PBS and mode-filtered through a Ib-ym pinhole to elim¬ 
inate the reflected light from the other interfaces, ex¬ 
cept for that of the trapped particle. We maximized the 
guided power by adjusting the trap beam alignment and 
the correction collar of the objective lens m- 

Figures 1(b) and 1(c) show the numerically calculated 
intensity distribution based on vectorial diffraction the¬ 
ory considering the following experimental parameters 
(as was discussed in our previous study m and refer¬ 
ences are therein): NA of the objective lens is 1.2, focal 
length f = 3 mm, filling factor = 1.0, vacuum wavelength 
Ao = 780 nm, and the refractive index (rim) of D 2 O is 
1.324. The intensity distribution is markedly different de¬ 
pending on the input beam polarization states: an elon¬ 
gated distribution for linearly (x) polarized beam (Fig. 
1(b)) and an isotropic distribution for circularly polar¬ 
ized beam (Fig. 1(c)). The minimum spot size (wq) was 
estimated to be 293 nm through a Gaussian fit of the 
numerically calculated field pattern nsun], whereas the 
conventional definition of optical resolution (the radius 
of the primary Airy disk) gives the relevant value as 270 
nm. Based on our experimental results, the beam waist 
(k/q) is chosen to be 300 nm, which is 2% larger than 
the ideal case, although we elaborated on producing the 
tightest spot size as possible. 

Part of the SV^ (red points) and SVy (blue points) data 
are presented in Fig. 2(a), and the corresponding power 
spectral density (PSD) is shown in Fig. 2(b) with the 
linearly (x) polarized trap beam. The beam power in 
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FIG. 2: (a) Part of QPD outputs: SVx (red) and SVy (blue), 
(b) PSD of QPD outputs: Q (red) and □ (blue) represent 
the PSD of SVx and SVy, respectively. Black and gray lines 
represent the noise levels of OT. The two vertical lines indicate 
the corner frequencies: fc,x (red solid line) and fc,y (blue 
dotted line). 
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the trap region was estimated to be Pq = 40 mW. The 
stochastic motion of the trapped particle can be mod¬ 
eled in the form of the Langevin equation with an addi¬ 
tional Hookean force term, where the restoring force in 
a specific direction is characterized by the correspond¬ 
ing stiffness ki^ where i G {x^y^z). The power spectra 
for each degree of freedom can be approximated to the 
Lorentzian function Si{f) = + /^) [ 20 ] . 

where D is the diffusion constant. The corner frequency 
/c,i is given by kij27 t^ and 7 is the frictional constant 
given by Stokes’ law, 7 = dirya. The dynamic viscosity 
(77) of D 2 O is 1.12 xlO ^ kg/m/s and the nominal radius 
of the PS spheres is a = 150 nm. The thick solid lines (red 
and blue curves in Fig. 2(b)) represent the results of the 
Lorentzian fit following the procedures reported in Ref. 
Enmi], which give the corner frequencies for each degrees 
of freedom in the transverse plane; {fc,x^ fc,y) = (2497 ± 
25, 3616 ± 33) Hz, respectively. The associated trap stiff¬ 
nesses {kx, ky) are (49, 71) pN//im and the corresponding 
stiffness asymmetry factor, st = 1 — kx/ky^ in the trans¬ 
verse plane is = 0.31 for the linearly {x) polar¬ 

ized trap beam. Our previous numerical study based on 
the revised generalized Lorenz-Mie theory (GLMT) pre¬ 
dicted the transverse asymmetry factor = 0.3 

under an aberrations-free condition [16], which is in a 
reasonable agreement with the experimentally measured 
value. 


III. CHARACTERIZATION AND CONTROL OF 
STIFFNESS ASYMMETRY 

The landscape of the optical potential depends not 
only on the EM filed distribution in the trapping re¬ 
gion, but also on the characteristics of the particle be¬ 
ing trapped, i.e., its size, shape, and index of refraction, 
and the index of refraction of the suspension medium 
as was demonstrated in the previous studies HDHUHS]. 
We focused on the systematic investigation of the polar¬ 
ization effect on the trap potential asymmetry using the 
polystyrene (PS) spheres with the diameters of 300 ± 
15 nm (< Ao) and 1000 ± 30 nm (> Aq). Polarization- 
induced stiffness asymmetry was investigated as a func¬ 
tion of the prepared input polarization state. The exper¬ 
imental parameters used in the preliminary experiment 
are used here again. 

The angular variation of the trap stiffnesses kx and ky 
are estimated from PSD measurements for the linearly 
polarized input beam rotated by Or using HWP2, where 
kx and ky are stiffnesses along the horizontal {x) and ver¬ 
tical (y) directions in the laboratory frame, respectively 
(referred to as the QPD segmentation axes). GLMT 
calculation results are compared with the experimental 
results in terms of the normalized stiffness, ki = ki/k^ 
where the average stiffness is defined as k = (kx ^ky)/2. 
For the 300-nm PS spheres presented in Fig. 3, the sym¬ 
bols (•, ■) and the gray (solid, dashed) curves denote 
the experimental and theoretical values of the normal- 



FIG. 3: Angular dependence of the normalized trap stiffness, 
kx and ky. Results of the 300-nm PS sphere: The filled gray 
(circle, square) symbols and gray (solid, dashed) curves rep¬ 
resent the experimental and theoretical values of kx and ky, 
respectively. Results of the 1-ym PS sphere: The empty red 
(circle, square) symbols and red (solid, dashed) curves rep¬ 
resent the experimental and theoretical values of kx and ky, 
respectively. 


ized transverse stiffnesses, kx and ky, respectively. The 
orthogonal pair of kx and ky shows periodic variation as 
the input beam polarization angle Or is rotated, which 
could be inferred from the rotation of the elongated in¬ 
tensity distribution in Fig. 1(b). The average value of 
the experimentally measured trap stiffness was ^exp = 
59 pN/ ym, and the GLMT calculation with the minimum 
spot size {n/Q = 300 nm) predicts ^glmt = 60 pN//im, 
showing good agreement with the experimental results. 
From the angular measurements, the stiffness asymme¬ 
try is estimated to be st,exp = 0.29, and the GLMT 
calculation gives st,glmt = 0.28. 

Figure 3 also contains the results of the 1-ym PS 
spheres: the symbols (O, □) and the red (solid, dashed) 
curves denote the experimental and theoretical values of 
the normalized transverse stiffnesses, kx and ky, respec¬ 
tively. We note the inversion of the angular stiffnesses, 
kx and ky, between the 300-nm and 1-ym PS spheres; 
for the 300-nm PS spheres, kx < ky at Or = 0 ° owing 
to the elongated intensity distribution of the linearly po¬ 
larized beam, whereas kx > ky in case of the 1-ym PS 
spheres with the same polarization configuration. This 
implies significant EM field redistribution by the pres¬ 
ence of a large dielectric particle, even in spherical shape. 
The average value of the experimentally measured trap 
stiffness and the theoretical estimation of the 1-ym PS 
spheres were (^exp, ^glmt) = (97, 125) pN//im, and 
the stiffness asymmetry factors were (st,exp, «§t,glmt) 
= (-0.1, - 0 . 12 ), respectively. Although the experimen¬ 
tal results and theoretical estimations show qualitative 
agreement (sign change of asymmetry factor), there is 
about 20 % discrepancy in the 1-ym PS spheres results, 
which probably originates from the fifth-order Gaussian 
beam method used in GLMT calculation [161 Ell l23) . 
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High-order Gaussian beam method provides the focused 
EM field distribution in the transverse plane with errors 
of within a few percentages; however, it renders a steeper 
intensity variation along the axial direction compared to 
the vectorial diffraction theory [24] . 

Figure 4(a) shows a demonstration of trap asymme¬ 
try control, which proceeds by adjusting the polariza¬ 
tion state of the input beam. The polarization state of 
the input beam is prepared as = Eo(cos6>qwp^ + 
i sin^Qwp^), where ^qwp is the angle between the slow 
axis of the QWP and the horizontal {x) axis. As the po¬ 
larization state becomes circular polarization at 6 >qwp = 
45°, the strong asymmetry of trap stiffness diminishes 
for both the 300-nm and l-/im PS spheres. In partic¬ 
ular, the inverted stiffness asymmetry (st < 0 ) of the 
l-/im PS spheres could also be readjusted to be isotropic 
with a circularly polarized trap beam, which can be in¬ 
ferred from the rotational symmetry of the isotropic field 
distribution of a circularly polarized beam and the PS 
spheres. To explore the trap potential landscape in the 
transverse plane, we investigated the angular variation 
in the balanced trap stiffnesses as follows. For each an¬ 
gles of HWP2, the direction of the input beam polar¬ 
ization is rotated to 6 >r = 26>hwp 2, and the QWP is 
rotated to produce a circularly polarized beam for the 
given polarization angle: ^qwp = + 45°. Figure 4(b) 

shows moderate variation in the normalized stiffnesses 
compared to the previously discussed angular variation in 
the trap stiffnesses for linearly polarized light (presented 
as light curves for comparison). The average stiffness 
for the 300-nm (l-/im) PS spheres was 59 (98) pN//im 
with the Icr-level standard deviation of 0.9 (0.9) pN//im, 
corresponding to fluctuations of 1.4 (1) %, respectively. 

Figure 5 shows the linear dependence of trap stiffness 
on the beam power for linear and circular polarization at 
a fixed trap depth of 50 /im from the inner wall of the 
cover glass. Other measurements (not presented) of trap 
stiffness as a function of the trap depth over the range 
of 10-150 /im show a few percentages of random fluctu¬ 
ations, as is typical for a water-immersion objective lens 
nans]. The trap stiffness per beam power for each de¬ 
gree of freedom is found by the linear fit of the experimen¬ 
tal data: = (1.21, 1.74, 0.27) pN//im/mW 

and (k'y.,k'y,k'^)^^ = (1.47, 1.48, 0.27) pN//im/mW for 
the 300-nm PS spheres; and {k'^.k'y.k'^)^^ = (2.58, 2.30, 
0.84) pN//im/mW and {k'y.,k'y,k'^)^^ = (2.47, 2.46, 0.83) 
pN//im/mW for the l-jam PS spheres, where the super¬ 
script denotes the polarization states. 


IV. CONCLUDING REMARKS 

The anisotropic field distribution of a highly focused 
EM fields manifests as dramatically modified optical po¬ 
tential landscape depending on the participating dielec¬ 
tric spheres, i.e., in terms of the sign and magnitude of 
stiffness asymmetry factor: = (0.29, -0.1) for 300- 
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FIG. 4: Trap stiffness control via the polarization state of 
trap beam. Results of the 300-nm PS sphere: The filled gray 
(circle, square) symbols and gray (solid, dashed) curves rep¬ 
resent the experimental and theoretical values of kx and ky 
respectively. Results of the l-/im PS sphere: The empty red 
(circle, square) symbols and red (solid, dashed) curves rep¬ 
resent the experimental and theoretical values of kx and ky, 
respectively, (b) Trap stiffnesses, kx and ky, variation of the 
balanced potential. Same symbols are used. The horizontal 
gray (red) lines represent the 1 cr-level standard deviation of 
the 300-nm (l-/im) PS spheres, respectively. 



(a) (b) 


FIG. 5: Trap stiffness as a function of trap beam power: (a) 
300-nm and (b) l-/im PS spheres. The (square, circle, trian¬ 
gle) symbols denote {kx,ky,kz), respectively, and the (filled, 
empty) status of the symbols distinguishes the polarization 
state (LP, CP), respectively. For most data points, the error 
bars (±1(7 level) are smaller than the size of the symbols. 
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nm and l-/im PS spheres, respectively. In case of the ad¬ 
justed optical trap realized by a simple polarization con¬ 
trol, the angular measurements of the orthogonal stiffness 
pair show less than 1.5% variation in a statistical sense 
(kr^xpl ^ 0.03). Our research translates the polariza¬ 
tion state of the trap laser beam, which is responsible for 
trap potential asymmetry, into a control parameter for 
tailoring the optical potential landscape of OTs in the 
transverse directions. 
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